K 1-6: an asymmetric planetary nebula with a binary central star 
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Abstract: 

We present new imaging data and archival multiwavelength observations of the little studied emission 
nebula K 1-6 and its central star. Narrow-band images in Ha {+ [N ll]) and [O ill] taken with the Faulkes 
Telescope North reveal a stratified, asymmetric, elliptical nebula surrounding a central star which has 
the colours of a late G or early K-type subgiant or giant. GALEX ultraviolet images reveal a very hot 
subdwarf or white dwarf coincident in position with this star. The cooler, optically dominant star is 
strongly variable with a period of 21.312 ± 0.008 days, and is possibly a high amplitude member of the 
RS CVn class, although an FK Com classification is also possible. Archival ROSAT data provide good 
evidence that the cool star has an active corona. We conclude that K 1-6 is most likely an old bona fide 
planetary nebula at a distance of ~1.0kpc, interacting with the interstellar medium, and containing 
a binary or ternary central star. The observations and data analyses reported in this paper were 
conducted in conjunction with Year 11 high school students as part of an Australian Research Council 
Linkage Grant science education project, denoted Space To Grow, conducted jointly by professional 
astronomers, educational researchers, teachers, and high-school students. 

Keywords: stars: AGB and post-AGB — biriarie: general — planetary nebulae: general — planetary 
nebulae: individual: K 1-6 



1 Introduction 

Planetary nebulae (PN) are an important, albeit brief 
(<10^ yr) phase in the evolution of low- to intermediate- 
mass stars (~1 to 8 M0). They are formed at the end 
of the asymptotic giant branch (AGB) phase, as the 
red giant ejects its outer envelope in a slow but massive 
'superwind'. After the envelope ejection, the remnant 
core of the star increases in temperature before nuclear 
burning ceases and the star fades to become a white 
dwarf. The hot surface temperature of the central star 
causes the previously ejected material to be ionized, 
which becomes visible as a PN, its shape influenced by 
the interaction between the old red giant envelope and 
a tenuous, fast wind from the hot central star (Kwok, 
Purton & Fitzgerald 1978). 

Recent work has shown a remarkable diversity in 
PN morphologies and central star properties (Frew & 
Parker 2010) . Since evidence is increasing that the 'PN 
phenomenon' is heterogeneous, and PNe are likely to 



be formed from multiple evolutionary scenarios, it is 
important to ascertain the nature of any poorly studied 
emission nebula that has been ambiguously classified 
in the literature. 

The object known as K 1-6 is a case in point. It 
was discovered by Kohoutek (1962), who noted a faint 
nebulosity visible on Palomar Observatory Sky Survey 
(POSS-1) red and blue prints. He classed it as a prob- 
able PN, even though he had no spectroscopic data, 
nor did he find a candidate ionizing star. Due to its 
low surface brightness, it has been little studied over 
the half century since its discovery. Perek & Kohoutek 
(1967) included it in their Catalogue of Galactic Plan- 
etary Nebulae under the designation PK 107-1-21.1. 
Cahn & Kaler (1971), assuming it was a PN, estimated 
a red flux from the POSS. Kwitter, Jacoby & Lydon 
(1988) took CCD images centred on the nebula, but 
also failed to find an ionizing star. Consequently, its 
status as a bona fide PN was in doubt. Acker et al. 
(1992) classed it only as a possible PN, with the com- 
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ment that it was a faint, highly reddened object show- 
ing weak Hq emission, based on unpubhshed spectra. 
On the other hand, Kohoutek (2001) included it as a 
bona fide PN, without elaborating on any reasons for 
its classification. 

K 1-6 is visible on both blue and red POSS-II 
plates, but reflection nebulae are also known to occur 
at high Galactic latitudes (e.g. Goerigk et al. 1983). 
Since it is a high latitude object at a very northerly 
declination (+74°) we cannot use the IPHAS (Drew 
et al. 2005) or SHS (Parker et al. 2005) surveys to 
investigate its morphology on deep Ha images. How- 
ever, it is visible on a low-resolution (96"), continuum- 
subtracted Ha image downloaded from the Virginia 
Tech Spectral line Surve£| (VTSS; Dennison, bimon- 
etti & Topasna 1998), from which we can conclude it is 
a genuine emission nebula. Since the latitude is quite 
high, an informed guess would class it as a PN rather 
than a conventional Hll region. However, deep, higher 
resolution, narrowband images are needed to better 
understand its morphology, as well as evidence that a 
low-luminosity ionizing star is present in the vicinity 
before it can be definitively classified as a PN. 

In this paper, we describe our efforts to ascertain 
the nature of K 1-6, selected as an object worthy of 
study as part of our Space to Grow project. This is 
a science education project conducted jointly by pro- 
fessional astronomers, educational researchers, teach- 
ers, and high-school students supported by an Aus- 
tralian Research Council (ARC) Linkage Grant. In 
Section [2] we give some background on this science 
education project, in Section |3] we describe new nar- 
rowband imaging data acquired with the Faulkes Tele- 
scope North, as well as archival multiwavelength data 
obtained from the Virtual Observatory (VO). We dis- 
cuss our interpretations on the nature of K 1-6 in Sec- 
tion 131 offer suggestions for future work in Section [S] 
before summarising our conclusions in Section [6] 

2 The Faulkes Telescope ARC- 
linkage Project 

The Space to Grow Faulkes Telescope Project is a 
three year ARC Linkage Grant that uses astronomy 
as the vehicle to engage students in real science dur- 
ing their final years of secondary school (Years 10- 
12). Launched in 2009, this major education project 
is administered by Macquarie University and follows 
the success of our Australian School Innovation in Sci- 
ence, Technology and Mathematics (ASISTM) pilot 
projec10 to help address the serious problem of poor 
student engagement and retention in science in Aus- 
tralian high schools. The industry partners in this 
ARC linkage project are Charles Sturt University, the 
Catholic Schools Offices of Parramatta and Bathurst, 
the NSW Department of Education (Western Region) 
and the Las Cumbres Observatory Global Telescope 
network (LCOGT)0 

^http: / /www. phys.vt.edu/~halpha/ 
^http: / /www. astronomy.mq.edu.au/asistm/ 
^http://Icogt.net/ 



The project is funded until 2012, and we will poten- 
tially have up to 40 schools, 200 teachers, and approxi- 
mately 5000 students from across different educational 
sectors in NSW involved in the project. Students and 
their teachers have the opportunity to work with pro- 
fessional astronomers on existing projects and create 
their own research proposal and submit requests to 
obtain data from the two, 2-metre Faulkes Telescopes 
located in Coonabarabran, NSW and Maui, Hawaii. 
These are the world's largest optical telescopes regu- 
larly available for science education projects in schools. 

A suite of detailed learning materials (accessible 
at the Space To Grow Websitcl) is being developed 
to guide students' learning and help them develop the 
necessary skills to both apply for telescope time and 
analyse the scientific data obtained from such obser- 
vations. In addition, ongoing teacher professional de- 
velopment is embedded within the project to support 
teachers in competently covering the astronomy con- 
tent of the NSW science curriculum and to give them 
the confldence to use technology with their students in 
science lessons. 

A similar project exists in the context of radio as- 
tronomy, the CSIRO's PULSEQParkes ProjeclQ where 
students have the opportunity to control the 64-m Parkes 
radio telescope to observe galactic pulsars, analyse their 
own data on these objects, and contribute to profes- 
sional astronomers' science projects (Hollow et al. 2008; 
Hobbs et al. 2009). The overall aim of these projects 
is to engage students in real science and to support 
teachers in dealing with the processes and approaches 
that motivate science students. Both the Space to 
Grow and the PULSEQParkes projects have embed- 
ded educational research procedures to investigate the 
impact of these science experiences on both students 
and teachers. 

The observations and data analyses reported in 
this paper were conducted in conjunction with Year 11 
high school students with the goal of involving them 
in real science. One of us (J.S.) is a science teacher 
who was originally involved in the ASISTM project 
that provided students with access to the Faulkes Tele- 
scopes. We emphasise that the students were embed- 
ded in the "journey of scientific discovery" throughout 
the study of this interesting nebula, with in-depth dis- 
cussion sessions conducted as part of school visits by 
two team members (D.J.F. and M.F.). More specifi- 
cally, the students generated the colour composite im- 
age of the PN using various software packages, and, in 
conjunction with their teacher (JS), used online tools 
to retrieve the various archival data sets utilised in the 
paper, before generating the light curves for the cen- 
tral star, estimating the period and amplitude, and 
investigating the variability type based on the avail- 
able data. We emphasise that all software packages 
used in the analysis are freely available to high school 
teachers and their students. 



*http: / /www. astronomy.mq.edu.au/space2grow/ 

^http: / /outreach. atnf.csiro.au/education/pulseatparkes/ 
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3 Observations 

3.1 Faulkes Telescope Images 

Narrowband images of the nebula were obtained with 
the 2-metre f/10 Faulkes Telescope North (FTN) at 
Haleakala Observatory, Maui. The Merope camera 
was used which contains an E2V CCD with 2048 x 2048 
pixels. The plate scale of 0.14 arcsec pixel"^ leads to 
a field-of-view of 4.7arcmin on a side. Two prelimi- 
nary 600 sec exposures in each of two narrowband fil- 
ters centred on the Ha A6563A and [O ill] A5007A lines 
were acquired on the night of 30 November 2009, with 
the goal of better understanding the nebular morphol- 
ogy. Note that the Ha filter (Aeff = 6560A, FWHM 
— 92A) transmits the two red [Nil] lines in the wings 
of the filter. Shorter exposures were also taken in the 
standard Johnson B and V filters. Prior to download- 
ing, the raw images were pipeline-processed at LCOGT 
which included standard flat-fielding and bias subtrac- 
tion routines (Lewis et al. 2010). The resultant images 
were then retrieved as 2x2 binned FITS format images 
from the LCOGT website. 

The observing conditions on the night were not 
optimal however, as the frames were taken through a 
large airmass with poor seeing, and due to the long ex- 
posure times at large zenith angle, the images suffered 
from modest smearing due to tracking error and/or 
wind. Consequently, we obtained an additional 27 five 
minute exposures in [Olll] on 17/18 April 2010 and 
15 five minute Ha exposures on 8 May 2010. Images 
were rejected if they had visible cosmetic defects or 
fiat-fielding problems. We finally used 18 individual 
[O III] frames and 10 Ha frames for generating deeper 
images with higher S/N. The individual images were 
aligned in the Astrometric43 software package using 
an object detection and matching algorithm, and then 
median stacked to eliminate bad pixels and cosmic ray 
hits. The resulting final images are presented and dis- 
cussed in i|4l 

3.2 Archival Data 

To supplement our imaging data, we searched for multi- 
wavelength archival data of the nebula and central star 
using Virtual ObservatorjQ tools, accessible to profes- 
sional and amateur astronomers, educators, and stu- 
dents. Most of the data were retrieved and visualised 
primarily via the Aladin Sky Atlas from the Centre 
de Donnees astronomiques de Strasbourg (CDS), and 
the SkyView Virtual ObservatorjQ. Oostra (2006) has 
discussed the utility of using online databases for as- 
tronomical research in the classroom. 

We used archival images from the Galaxy Evo- 
lution Explorer (GALEX) satellite to look for a hot 
source in the vicinity, as our central star candidate is 
too cool to ionize the nebula. GALEX is a space-based 
50 cm telescope currently undertaking a nearly all-sky 
imaging survey using two ultraviolet (UV) filters: an 
FUV band (covering 1350-1750 A) and an NUV band 

'^http:/ /www.astrometrica.at/ 

'^http:/ /www. ivoa.net/ 

*http: / /skyview. gsfc.nasa.gov/ 



(1750-2750 A). Further details of the survey are given 
in Martin et al. (2005) and Morrissey et al. (2007). 
We visualised images downloaded via the GalexView 
1.4 web-based utilitjQ. A prominent, very hot star 
(FUV flux > NUV flux) was noted on Field AIS.46, 
designated GALEX J200414.5-F 742535 (see Figure [l|. 
The coordinates of this source from GALEX are within 
2" of our candidate central star. The positions are 
formally identical as the GALEX angular resolution is 
4.5-6" (Martin et al. 2005). We emphasise that no 
other conspicuous FUV source is visible inside a 15' x 
15' region centred on K 1-6 (see Figure [l]), so this star 
is clearly the ionizing source for the emission nebula. 
Our central star candidate is the optically brightest 
star visible in the interior of the emission nebula seen 
on our FTN images, and as we show below, has the 
colours of a G- or K-type star. We adopt its posi- 
tion from the Naval Observatory Merged Astrometric 
Dataset (NOMAD j3 accessed through the VizieEQ 
service at the CDS. We also adopt these coordinates 
as the nebular position (see Table [T}. 

Once the central star was identified, additional data 
was retrieved using VizieR, using a 5" search radius 
around the nominal position from NOMAD. Additional 
photometric data was obtained from the Tycho-2 cat- 
alogue (H0g et al. 2000) and the Two Micron All 
Sky Survey (2MASS; Cutri et al. 2003; Skrutskie 
et al. 2006), supplemented by lower-precision magni- 
tudes from the various scanned Schmidt plate surveys. 
Time-series photometric data were obtained from the 
Northern Sky Variability Survey websitfQ (NSVS; Wozniak 
et al. 2004) and The Amateur Sky Survey (TASS; 
Droege et al. 2006) webpageQ- 

We also noted an X-ray detection in the ROSAT 
All-Sky Bright Source Catalogue (Voges et al. 1999) 
with the designation IRXS J200413.0-I- 742533 (FigureH]). 
This object is coincident in position with both the 
cool star seen on the DSS and the GALEX source to 
within the nominal errors. The source was observed by 
ROSAT PSPC for 2350 sec (Seq. ID RS930419N00) at 
a count rate of 0.035 ± 0.006 counts s~^, thus resulting 
in 83 ± 15 backgound-subtracted counts. Further im- 
plications of the associated X-ray source are discussed 
in Section [4.51 

Finally, we note the absence of a nebular radio de- 
tection in the survey of Zijlstra, Pottasch & Bignell 
(1989), who gave an upper limit of 4mJy at 1.5 GHz. 
We visualised a 20' NRAO VLA Sky Survey (NVSS) 
image (Condon et al. 1998) using SkyView and con- 
firm this result. An unrelated background source is 
located 2.3' north-west of the nebular position. 



^http: / /galex. stsci.edu/galex View/ 
http: / / www. nofs . navy, mil/ nomad 
-'^-'^http: / /vizier, u-strasbg.fr/cgi-bin/ VizieR 
-'-^http://skydot.lanl. gov /nsvs/nsvs.php 
-"^^http: / /www. tass-survey.org/ 
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Figure 1: Left: Red DSS2 image of K 1-6, showing the faint nebula at centre. Middle: Corresponding 
GALEX NUV image of the same field. Right: Corresponding GALEX FUV image, showing the prominent, 
very hot ionizing source. Each image is 15' on a side, with north at top. 



4 Analysis 

4.1 Nebular Morphology 

We present our narrowband and colour-composite im- 
ages of K 1-6 in Figure [2] The Ha + [N ll] image shows 
the nebula to have an asymmetric appearance, with a 
brighter rim on the south side, while the [O ill] image 
shows the O^^ emitting zone to be roughly circular 
and restricted to a zone around the best candidate 
for the central star, i.e. the nebula shows ionization 
stratification. Similar stratification is commonly seen 
in evolved PNe with faint central stars on the white 
dwarf cooling track. The dimensions of the nebula 
were measured from the combined FTN Ha image as 
198 X 160". The available observational data on the 
nebula are summarised in Table [T] 

The position of our candidate central star is offset 
from the nebular centroid, closer to the brighter rim of 
the nebula as seen in Ha -I- [Nil] light (see Figure [21). 
We explain this arcuate morphology and the location 
of the ionizing source as a product of an nebular inter- 
action with the interstellar medium (ISM) as discussed 
by Wareing et al. (2007), Wareing (2010) and Sabin 
et al. (2010). We discuss the PN/ISM interaction in 
more detail in Section [4.61 

On our short-exposure FTN B images, the cool 
star appears perfectly round, showing that the hot star 
must have a separation of <1.2" (if it has V~15.8; see 
below). There are also five unrelated field stars within 
30" of the central binary; the closest has V ^16 and 
is 7.2" south-east (see Figure [2]). Since the positions 
of the cool and hot components of the central star are 
identical within the errors, we can state that the op- 
tically cool star has a very strong UV-excess and is 
potentially part of an unresolved binary system. In 
addition, we find the cool component to have a sig- 
nificant proper motion. The values from the USNO 
CCD Astrograph Catalog (UCAC3; Zacharias et al. 
2010), the Tycho-2 catalogue and the USNO-B cata- 
logue (Monet et al. 2003) are highly consistent, espe- 
cially in the S component. The values are much larger 



than the formal errors, and we take the detection to 
be real and not spurious. We adopt the UCAC3 value 
oi fj, = 29.1 masyr"^ in position angle 179°. 

This motion is practically due south, and is per- 
fectly consistent with the morphology of the bow shock 
on the southern edge of the nebula seen on our FTN 
Ha images. This is a key piece of evidence to support 
our hypothesis that the cool star is a physical com- 
panion to the true ionizing star. While the GALEX 
astrometry is not precise enough to rule out a chance 
superposition, we consider it improbable for the cool 
star to have a proper motion vector in the correct ori- 
entation, without it being associated with the PN. The 
asymmetry in the PN shell also leads to the conclusion 
that the nebula was ejected from the central binary 
system, and is not ionized ambient ISM (see Frew & 
Parker 2010). We assume hereafter that this is an un- 
resolved binary system, where the hot component is 
ionizing the surrounding nebula. Recall that we place 
an upper limit of 1.2" on the separation of the compo- 
nents from our FTN B short exposure images. 

4.2 K 1-6, a bona fide PN 

Is K 1-6 is a highly evolved PN? We note that PN- 
like nebulae with parabolic bow-shock morphologies 
are known (see Frew & Parker 2010, for a review). 
These include pulsar-wind nebulae and bowshocks as- 
sociated with nova-like cataclysmic variables (CVs), 
e.g. EGB 4 (Ellis, Grayson & Bond 1984; Krautter, 
Klaas & Radons 1987; HoUis et al. 1992) and Fr 2- 
11 (Frew, Madsen & Parker 2006; Frew et al. 2010b). 
The peculiar nebula Abell 35 (Jacoby 1981) has a very 
similar morphology to Fr 2-11 and is now thought to 
be a bowshock nebula inside a photoionized ambient 
H II region (Frew 2008). In both these latter cases, the 
bowshock is only visible in the light of [O ill] and is 
interior to the external edge of the ionized region. The 
observed morphology of K 1-6 is different and simi- 
lar to other one-sided PN; recall Figure [2] which shows 
that K 1-6 has an exterior bowshock which implies the 
ionized nebula is moving with respect to the ISM (cf. 
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Figure 2: Left: FTN Ha image of the one-sided, evolved PN K 1-6. Middle: FTN [O ill] image of K 1-6. 
Right: Colour image combined from Ha (red), [Oiii] (blue) and combined Ha -I- [Oiii] (green). The 
central star is marked. The closest of several unrelated field stars is 7.2" SE of the central star, and partly 
blended with it on these images. Each image is 4.7' across, with north-east at top left. 



Table 1: Basic properties of the nebula, K 1-6. 



a (J2000) 


20 04 14.25 


6 (J2000) 


4-74 25 36.4 


I 


107.036 


b 


-1-21.384 


Dimensions 


198" X 160" 


log F(Ha) (cgs) 


-11.32 ± 0.07 


E{B - V) 


0.22 


ch/3 


0.32 


Distance 


~1.0 kpc 


2-height 


~365 pc 



the emission nebula ionized by PHL 932; Frew et al. 
2010a). 

We have calculated the ionized mass of the nebula 
using equation 7 from Hua & Kwok (1999), which uses 
the measured diameter and integrated Ha flux, and an 
estimate of the distance, taken from Section [4.41 

Mion = 0.032 {e/0.6f-''e^'^D^''Fo{Ua)°'^ Mq (1) 

where 9 is in arcmin, D is in kpc, and Jb(Ha) = 
F{Ua) X iQO-e^xcHp jg ^jji^g of iQ~i2 ergcm^^s-^ 

We assumed a volume filling factor, e = 0.6 follow- 
ing Hua & Kwok (1999), and determine the integrated 
Ha flux in the following manner. No PN flux stan- 
dards were observed, so we perform aperture photom- 
etry from the calibrated VTSS Ha image (Dennison et 
al. 1998), following Frew, Parker & Russeil (2006). We 
measure a Ha flux of 4.8 ± 0.8 x 10~^^ ergcm~^s~^ 
through a 7' aperture. No correction for [Nil] emis- 
sion is needed as these lines are not passed by the nar- 
rowband VTSS filter, but we have corrected for the 
known zero-point error of the VTSS images (Gaustad 
et al. 2001; Finkbeiner 2003). The log fiux value is 
presented in Table [1] 

At a distance of 1.0 kpc, the ionized mass of K 1-6 
is 0.15 Mq, consistent with other PNe, which typically 



have ionized masses ranging from 5x10"^ to 3 Mq 
(Frew & Parker 2010). It is far too massive to be 
an evolved nova shell, which usually has a mass of 
the order of 10"* Mq (Cohen & Rosenthal 1983). In 
addition, the observed morphology is quite unlike that 
of a nova shell. 

Based on our new images, K 1-6 shows no evi- 
dence of any bipolarity or axisymmetric, low-ionization 
structures (LIS; Gongalves, Corradi & Mampaso 2001) 
which may be indicative of a close binary nucleus (see 
Miszalski et al. 2009, 2010a). The infiuence of bi- 
narity on PN formation and shaping is an important 
topic in the field, and is currently being debated (see 
De Marco 2009 for a recent review). Hence, any new 
binary central stars with periods in the weeks to years 
range are very important. Since we do not yet have 
any spectroscopic data on the central star, we have no 
information on the separation and orbital period of the 
binary components, but it is tempting to identify the 
21.3 day photometric period (see Section |4]3TT} as the 
orbital period of the system. However, such a system 
would be unusual in a non-bipolar PN. 

In summary, our FTN narrowband imaging data 
reveal a faint elliptical nebula that shows an asym- 
metric morphology indicative of an interaction with 
the ISM. The available data suggest it is best inter- 
preted as an evolved, low-surface-brightness PN based 
on its bow-shocked morphology, ionization structure, 
and ionized mass. Following the nomenclature intro- 
duced by Acker et al. (1992), the nebula is designated 
as PN G107.0-f21.3. 

4.3 Stellar photometry 

We summarise the available literature photometry of 
the binary central star in Tabled In order to correct 
the observed colours of the cool star for reddening, a 
colour excess of E{B — V)— 0.22 mag is adopted from 
the dust maps of Schlegel, Finkbeiner & Davis (1998) 
on this sightline, though this is technically an upper 
limit. The visual absorption is then A„ = 0.68 mag 
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Table 2: Archival multi-wavelength photometry of 
the central star of K 1-6. 



Band 


Magnitude 


Source 


IP 


1 Q Q7 _|_ n ni 


OAT ,RX 


ivuv 


1 4 Ko -1- n ni 

i'i.UiJ _1_ U.W-L 




Bjp 


13.00 ± 0.29 


Tycho-2 


Vt 


12.12 ± 0.19 


Tycho-2 


V 


12.16 ± 0.19 


GSC2.3 


V 


12.48 ± 0.36 


TASS-IV 


R 


12.29 ± 0.22 


NSVS 


R 


12.15 ± 0.08 


UCAC3 


Ic 


11.23 ± 0.20 


TASS-IV 


J 


10.619 lb 0.021 


2MASS 


H 


10.009 ± 0.018 


2MASS 


Ks 


9.809 ± 0.015 


2MASS 


Ftjv — NuY 


-0.65 


GALEX 


Bt - Vt 


0.88 


Tycho-2 


V-Ic 


1.25 


TASS 


J - H 


0.61 ± 0.03 


2MASS 


H-Ks 


0.20 ± 0.03 


2MASS 



using the R — 3.1 reddening law of Cardelli, Clayton 
& Mathis (1989). 

The optical and 2MASS near-IR colours for the 
central star are quite red, consistent with a G/K giant, 
or an early K-type main sequence star. Since E{V — /) 
= 1.28 E{B - V) = 0.28, we can correct the mean V-I 
colour from TASS to {V-I)o = 0.97. This corresponds 
to a spectral type of Kl-2. Similarly, a reddening cor- 
rected J — H colour of 0.55 corresponds to a spectral 
type of K3-4. The Tycho data suggest an earlier type 
around mid-G, but the Br-magnitude might be influ- 
enced by the underlying hot continuum of the ionizing 
star. Based on all available data, we adopt a spec- 
tral type of G8-K3 for the cool star hereafter. The 
luminosity class is indeterminate until spectra are ob- 
tained. 

We also derived an estimate for the visual mag- 
nitude of the hot star visible on GALEX images in 
the following manner. Using the GALEX Fuv {^oH 
= 1528A) and Auv (Aeff = 2270A) magnitudes for 24 
known PN central stars that show little or no extinc- 
tion (mean T^f! = lOOkK), we estimate {F\jy — N\jy)o 
= -0.75 ± 0.14 and (Auv - V)o = -1-8 ± 0.2, after 
correcting the GALEX magnitudes with a standard ex- 
tinction law (Cardelli et al. 1989). We keep the _Fuv 
and iVuv magnitudes in the AB system. Our derived 
values are consistent with the synthetic colours given 
by Venues, Kawka & Nemeth (2010). 

Using the de-reddened Fw and Auv magnitudes 
for the central star of K 1-6, we estimate Vb — 15.1 and 
re-redden to get V ~ 15.8. This is ~3.3 magnitudes 
fainter than the cool star, which explains why there is 
no clear signature of the hot star in the optical and 
near-IR colours of the system. The resulting absolute 
magnitude assuming D = 1.0 kpc (see Section [4.4p is 
Mv — +5.1, consistent with the absolute magnitudes 
of central stars in other PNe (Frew & Parker 2010a). 
This magnitude places the ionizing star near the top of 
the white dwarf cooling track, consistent with a low- 
mass central star in a highly evolved PN. 
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Table 3: Estimated period and amplitude of vari- 
ability 



Sample 


n 


Range 




Period 


A 






MJD 




(day) 


(mag) 


NSVS 


604 


51274 - 


51633 


21.32 ± 0.05 


0.68 


TASS 


169 


52786 - 


54024 


21.31 ± 0.05 


0.72 


All 


773 


51274 - 


54024 


21.312 ± 0.008 


0.70 



4.3.1 Stellar variability 

Interestingly, the cool star is variable, based on data 
from the Northern Sky Variability Survey (NSVS; Woz- 
niak et al. 2004), where it is designated as NSVS 1225362, 
The variability was discovered by Usatov & Nosulchik 
(2008), who noted a range oi R = 11.97 to 12.74 with 
P = 21 days. These authors noted the associated X- 
ray detection and suggested the star was a likely RS 
CVn variable. 

We have independently determined the period us- 
ing both the available NSVS and TASS data. Ini- 
tially the data from NSVS was folded on an initial 
estimate of the period from the average separation be- 
tween apparent peaks in the lightcurve. New periods 
around this value were trialed and the amount of scat- 
ter around the mean curve were visually estimated. 
The point at which the scatter was minimised was de- 
termined to be the most likely period, 21.3 days. The 
error in this method was determined to be not less 
than 0.1 day due to its subjective nature. 

A more formal approach was then undertaken. We 
combined the TASS V and / data to form a synthetic 
red magnitude, which was found to be offset from the 
NSVS i?-band data by 0.40 mag, in the sense the TASS 
data are fainter. After correcting to a common zero 
point, the TASS and NSVS observations were com- 
bined to provide a data set of over 770 points span- 
ning 2749 days. We used the PeranscQ software pack- 
age using an ANOVA method to determine a period 
of 21.312 ± 0.008 days. The periods estimated from 
each data set are given in Table O 

The combined light curve as a function of phase is 
given in Figure [3l The phased- up light-curve is quite 
regular and is quasi-sinusoidal with an amplitude of 
~0.7 mag, but we note a slower rise to maximum fol- 
lowed by a faster decline. We also looked at the TASS 
V ~ I colour index as a function of phase. The star 
clearly gets redder as it fades (see the lower plot in 
Figure [3|, consistent with either a Cepheid or RS CVn 
interpretation. 

There is no evidence for any eclipses in the avail- 
able archival data. We also found no evidence for any 
rapid variations (on the order of minutes to tens of 
minutes) from differential photometry of the cool star 
on our FTN narrow-band and broadband B and V im- 
ages. In order to constrain the variability type of the 
cool star, we considered the following possibilities: 

• Irradiation variable - unlikely as amplitude seems 
^*http://www.peranso.com/ 



Frew et al. 



7 




0.0 0.5 1.0 1.6 2.0 



Phase 

Figure 3: Phased light curve of the central star of 
K 1-6 assuming a period of 21.312 days. The top 
graph shows all data points from NSVS (wideband 
R magnitude) and TASS (synthetic red magnitude 
— see text). The bottom curve shows the TASS 
V — I colour index as a function of phase. The star 
clearly becomes redder as it fades. 



to be too large for period (see De Marco, Hillwig 
& Smith 2008); 

• Ellipsoidal variable - unlikely as amplitude is 
too large, and light-curve is asymmetric; 

• Small-amplitude red variable (SARV) - these 
are red giants of types late K and early M; our 
amplitude and period are consistent but our star 
appears too hot. 

• Ordinary pulsating G or K subgiant or giant - 
rejected as the observed amplitude is much too 
large and the period too long; 

• Type II Cepheid - appears consistent with ob- 
served parameters, but calculated distance de- 
rived from the P-L relation is too large; 

• Anomalous Cepheid - these are metal poor stars 
and represent the extension of classical Cepheids 
to lower metallicities and masses - rejected as 
our period is too long; 

• Cataclysmic variable - rejected as there is no 
evidence of rapid variations on our FTN images, 
or any eruptive behaviour from archival data; 

• Classical symbiotic star - rejected as there is 
no evidence for high-density nebular gas (but 
we lack spectra) or eruptive behaviour; observed 
period appears too short. 

• Evolved, chromospherically active variable, e.g. 
RS CVn or FK Com type - also feasible but 
our star has an unusually large amplitude; hard 
X-ray emission is diagnostic. 

At first glance the amplitude, period and dered- 
dened V ~ I colour of the star appear to be consistent 
with it being a Type II Cepheid. The light curve has 
a slow rise and a faster decline which is unusual, but 



not unprecendented amongst this class (see Figure 3 of 
Soszyhski et al. 2008). The period falls in the distri- 
bution gap between the W Virginis and RV Tauri sub- 
classes, but the light curve is more like a W Vir star, 
as it lacks the alternating deep and shallow minima 
characteristic of the RV Tan class (Wallerstein 2002). 

However, the detection of hard X-rays (see sec- 
tion 14. 5p strongly suggests that it is a coronally ac- 
tive cool star, such as an RS CVn system, as first 
suggested by Usatov & Nosulchik (2008). It may in 
fact be a member of the V471 Tan subclass, or even 
an FK Com star, pending spectroscopic analysis. We 
note however that the observed amplitude (~0.7 mag) 
is quite large for either a FK Com or RS CVn star. 

4.4 Distance 

If the cool component seen in K 1-6 was a Type II 
Cepheid, the distance would be easy to estimate since 
these stars follow a well defined period-luminosity (PL) 
relation, which is 1.5-2.0 mags fainter at a given period 
than the classical Cepheid PL relation (Soszyhski et 
al. 2008). On the assumption that the star is a W Vir 
star, we calculate a distance using the mean reddening- 
corrected l/-band magnitude from the TASS and the 
relation of Pritzl et al. (2003). We obtain a distance of 
6.3 ± 0.6 kpc. However, we immediately discount this 
interpretation because the resulting stellar transverse 
velocity would be unphysical (~900 kms~^) and the as- 
sociated PN would be impossibly massive (>lOAf0). 
Furthermore, the z-distance from the Galactic plane 
would be 2.3 kpc, placing the PN well outside the dust 
and gas layer in the disk (Spitzer 1978; Dickey & Lock- 
man 1990), hence a strong ISM interaction would not 
be expected. In addition. Type II Cepheids are rare 
stars and the chance of having one in a physical binary 
with a PN central star is remote. Lastly, we note that 
Type II Cepheids are not strong X-ray sources. 

While we have discounted the cool star as a Type II 
Cepheid, we can still estimate an approximate photo- 
metric parallax from a combination of the reddening- 
corrected V — I and J ~ H colours. The Tycho-2 data 
are not reliable at this magnitude and we also decided 
against using the ifs-band magnitude from 2MASS, as 
it may be affected by hot dust — note that the star 
appears to have a weak detection in archival 12 /^m 
IRAS data, but it is not listed as an IRAS source. 

A preliminary estimate of the spectral type from 
the reddening-corrected colours is G8 - K3 (see above). 
However, we cannot derive the luminosity class from 
the available data. Since this is uncertain, a range 
of possible luminosity classes is given in Table 2] We 
adopt the absolute magnitudes for each class from Schmidt- 
Kaler (1982) and Cox (2000) in column 2. Using the 
reddening-corrected mean Vb = 11.8 ± 0.2, we de- 
rive the spectroscopic distances in column 3. and the 
height above the Galactic plane in column 4. We cal- 
culate the transverse velocity (in kms~^) for each case 
using 4.74 ^ D, where D is in pc and /i is in arcsec yr~^. 
These are given in column 5. We can immediately rule 
out the cool star being a main sequence star (the ve- 
locity is too low to form a pronounced bow shock in 
the PN) or a bright giant of luminosity class II (trans- 
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Table 4: Inferred properties for the cool compo- 
nent of K 1-6 as a function of spectral class. 



Spectral type 


Mv 


D 


z 


vt 






(pc) 


(pc) 


(kms-i) 


K3 V 


+6.7 


105 


38 


14 


Kl V 


+6.2 


130 


48 


18 


G8/K0 IV 


+3.1 


550 


200 


75 


KO III-IV 


+1.9 


960 


350 


130 


G8 III 


+0.8 


1580 


580 


220 


KO III 


+0.7 


1660 


605 


230 


K3 III 


+0.3 


2000 


730 


275 



verse velocity is unphysical, as for the Type II Cepheid 
interpretation). A luminosity class of either III or IV 
is consistent with the available data. 

We can use an alternative approach to find the dis- 
tance, based on the properties of the nebula, assuming 
it to be a bona fide planetary. Our recently developed 
Ha surface brightness-radius (SB-r) relation (Frew & 
Parker 2006; Frew 2008) uses the reddening-corrected 
Ha surface brightness to infer the radius of the PN (in 
pc). The angular size of the PN then directly gives the 
distance. Adopting the dimensions, Ha fiux, and red- 
dening of the PN from Table [T] we derive a distance of 
1.6 kpc using the preliminary SB-r relation of Pierce 
et al. (2004), or 1.1 kpc using a short trend applica- 
ble to lower- mass shells (Frew & Parker 2006, 2007). 
If we factor in representative 1-sigma uncertainties on 
the SB-r distances (Frew 2008), the extreme distance 
range is 0.8-2.0 kpc. Since the X-ray properties (see 
Section I4.5|l suggest the nearer end of this range is 
more likely, we adopt a working distance of 1.0 kpc for 
K 1-6 hereafter, and the inferred luminosity class for 
the cool star is then III-IV. 

4.5 X-ray analysis 

The central star of K 1-6 coincides with the ROSAT 
source IRXS J200413.0+742533 (Figure B]). We ex- 
tracted the ROSAT X-ray spectrum using the HEA- 
SOFT software package (Figure [5]). The background- 
subtracted X-ray spectrum shows a prominent peak at 
0.7-0.9 keV. Adopting a column density iVH = 1.3xlO^^ 
cm~^ consistent with the assumed E{B — V) = 0.22, 
the best-fit, thin plasma MEKAL model (Mewe, Kaas- 
tra & Liedahl 1995) has a temperature kT = O.Slg 2 
(or 6x 10^ K). We find that the soft X-ray component is 
poorly fitted at the adopted extinction, so the asymp- 
totic SFD reddening may be an overestimate. 

The unabsorbed X-ray fiux in the 0.3-2.0 keV en- 
ergy band is 3.1x10"^^ erg cm~^ s~^ with an un- 
certainty within a factor of two. At a distance of 
1.0 kpc, it implies an X-ray intrinsic luminosity of Lx 
= 3.7xl0^^(d/kpc)^ erg s"^ This luminosity is con- 
sistent with an RS CVn or FK Com star, but near 
the upper limit for these classes (Walter & Bowyer 
1981; Makarov 2003; Gondoin 2005). We also note 
that the luminosity is somewhat too bright for the 
measured kT after referring to the log Lx versus log T 
plot of Schmitt et al. (1990). This might also be an 
indication that the adopted distance is somewhat too 



large and/or the extinction has been overestimated. 

Assuming our X-ray temperature to be correct, we 
can also use it to place constraints on the nature of 
the binary star. The temperature appears too low to 
be a CV with an accretion disk (Eracleous, Halpern 
& Patterson 1991), and even though the combination 
of low kT and high Lx overlaps with symbiotic stars 
(e.g. Murset, Wolff & Jordan 1995), there is as yet 
no corroborating observational evidence that this is 
indeed a symbiotic system. The nebular morphology 
is also unlike the strongly axisymmetric (toroidal and 
bipolar) forms seen in true symbiotic outfiows (Corradi 
2003). However, spectroscopy of the central star might 
reveal high-density emission lines indicative of the new 
class of EGB 6-like nuclei (Frew & Parker 2010). 

The absolute bolometric magnitude of the cool com- 
ponent assuming D — 1.0 kpc is Mboi = +1.4 (equiva- 
lent to Lboi = 22L0). The distance- independent ratio, 
Rx = log (La;/Lboi) is ~ —3.4 which is close to the sat- 
urated level of X-ray emission from an active corona 
(Giidel, Guinan & Skinner 1997; Giidel 2004; Gtidel & 
Naze 2009). 

Our observed photometric modulation has a 21.3 
day period fSection I4.3.1|) . but it remains uncertain 
what this period represents. Assuming Tcft ~5000 K 
for the cool star, the observed luminosity leads to a 
radius of ~6.5 Rq. Furthermore, if we assume a mass 
of ~1.0-1.5 Mq, the fundamental pulsation period is 
calculated to be only 0.4-0.6 days, so pulsation does 
not appear to be the explanation of the observed light 
variation. 

Since the cool star appears to be a subgiant or 
giant, we cannot use a standard rotation-activity rela- 
tion applicable to main-sequence stars (e.g. Pallavicini 
et al. 1981; Giidel et al. 1997) to determine the rota- 
tion period from the Lx value. However, the observed 
period is consistent with the observed rotation periods 
of other coronally active cool giants. Gondoin (2005, 
2007) has shown there is a relationship beween the X- 
ray surface fiux and rotation for cool giants, both sin- 
gle and binary. Gondoin (2007) adopts the following 
relation: 

log(F,) = -0.95 log(P,ot) + 7.1 (2) 

For the cool star in K 1-6, the X-ray surface fiux is 
~10^ ergs~^cm~^, good to a factor of two. The pre- 
dicted rotation period using equation[2]is then only ~1 
day, which must be in error as this exceeds the critical 
rotation rate. Better data, including an independent 
determination of the star's rotation period, are needed 
to see if the star follows the relationship between activ- 
ity and rotation. We note that two active giants with 
comparable log -Fx values, FF Aqr and 39 Get, have 
rotation periods of 9.2 and 77.2 days respectively, and 
fall well away from the relation of Gondoin (2007); cu- 
riously, both have hot subluminous companions. 

In many respects this system resembles the old, 
high-latitude PN LoTr 5 (Longmore & Tritton 1980; 
Graham et al. 2004), which holds a binary nucleus 
with a hot WD and a chromospherically active, opti- 
cally variable G5III companion, IN Comae, which is 
a strong X-ray emitter (Guerrero et al. 2000; Montez 
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Figure 4: Smoothed ROSAT PSPC images of the field containing K 1-6 in the energy ranges 0.1-1.3 keV 
{left), 0.1-0.45 keV [center], and 0.55-1.3 keV [right). The position of K 1-6 is marked. The adjacent 
source towards the east is likely to be a background AGN, as suggested from its energy hardness. North 
is at the top, and east to the left. The scale bar on the first image represents 10'. 



et al. 2010). However, the optical amplitude of the 
variable cool component in K 1-6 is much larger. 




Figure 5: ROSAT PSPC spectrum and best-fit 
model for the central star of K 1-6. The spectrum 
can be deconvolved into two components: a softer 
component corresponding to the high energy tail 
of a hot stellar photosphere, and a harder compo- 
nent that is associated with emission from the cool 
active companion. 



4.6 PN/ISM interaction 

Wc note several PNc with similar stratified, one-sided 
morphologies to K 1-6. Examples of such 'sosio' PNc 
are Lo 10 (Schwarz, Corradi & Melnick 1992), BV 5-2 
(Manchado et al. 1996), Te 3 (Jacoby et al. 2010), 
Is We 1 (Xilouris et al. 1996; Tweedy & Kwitter 1996), 
and Abell 31 (Abell 1966; Tweedy & Kwitter 1994). 
The morphologies of these nebulae arc attributed to an 
ISM interaction in each case. We show now that this 
explanation is valid to explain the morphology seen in 
K 1-6. 

At our preferred distance of 1.0 kpc, the transverse 
velocity of the binary central star is 140 kms~^, which 
is at the high-velocity tail of the velocity distribution 

for the local thin disk population (Skuljan, Hcarnshaw 
& Cottrell 1999). The true space motion might be 
higher (the radial velocity is currently unknown) so it 
might also be a thick disk object or a rare halo PN. Ei- 
ther way, the high velocity explains the morphological 
asymmetry of the PN, due to a strong interaction with 
the ISM. Borkowski, Sarazin & Sokcr (1990) showed 
that significant distortion of a PN becomes apparent 
once the pressure of the expanding PN becomes equal 
with the ram pressure of the surrounding warm ISM. 
Following Borkowski et al. (1990) and Jacoby & Van 
de Steene (1995), this relation is expressed mathemat- 
ically as: 



no [Vt + Ve^pf = He 



(3) 



For our preferred distance of 1.0 kpc, the z-hcight 
is 365 pc. Using the best-fit distribution of HI with 
height in the Galaxy (Dickey & Lockman 1990), the 
ISM density is no = 0.06 cm ~^ at this ^-height. We 
have no kinematic information on K 1-6, so we as- 
sume an expansion velocity of 20 kms~^ for the shell 
(Weinberger 1989), and the root-mean-square electron 
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density of the PN is calculated as Ue — 20 cm ~^ based 
on the observed diameter and flux. 

Re-arranging equation [3] and solving for v, , we in- 
fer a space velocity of 350 kms~^ which is consistent 
to roughly a factor of two with the estimated trans- 
verse velocity at this distance. At smaller distances 
the ISM density is higher, and the required space ve- 
locity is reduced (but remember the transverse veloc- 
ity scales with D). Alternatively, the ISM is denser 
than predicted by Dickey & Lockman (1990) along this 
sight-line. Of course the central star could also have 
a significant radial velocity of a similar magnitude to 
the transverse component, in which case the star be- 
longs to the halo. Considering the assumptions of this 
rather coarse approach, we conclude that the observed 
velocity of the star is enough to explain the observed 
nebular interaction with the ISM at this z-height. 

5 Future Work 

Initially, medium-resolution, optical spectroscopy of 
the central binary (or ternary) star of K 1-6 is required 
in order to remove any ambiguity in its classification, 
and to determine its metallicity. Ultimately, an accu- 
rate spectral type and luminosity class for the central 
star will allow an improved photometric parallax for 
this interesting PN. If the distance is confirmed, then 
K 1-6 is an addition to the solar neighbourhood PN 
census (Frew & Parker 2010). This analysis might be 
hinting that there are other nearby, low-mass PN re- 
maining to be discovered, especially well off the Galac- 
tic plane. 

A combined time-series photometric and radial ve- 
locity study will put constraints on the physical pa- 
rameters of the central binary system. If the cool star 
is an RS CVn variable, we would predict strong Call 
H and K emission and Mgll emission which are sig- 
natures of strong chromospheric activity. Higher reso- 
lution spectroscopic observations would provide direct 
evidence for the cool star having experienced spin-up. 
Fast rotation in an old star could be a signature of 
mass accretion from a former companion in a detached 
binary system, now seen as the ionizing source of the 
PN. Jeffries & Stevens (1996) have proposed that a 
detached secondary in a wide binary can accrete part 
of the slow massive wind from an AGB companion. 
The gain in mass and angular momentum can spin up 
the secondary to a state of very rapid rotation; as the 
AGB star evolves to a white dwarf, the companion now 
appears as a wind-accretion induced rapid rotator (or 
WIRRing star). Soker & Kastner (2002) showed that 
this is viable on theoretical grounds, for cases where 
the initial separation was less than 30-60 AU. 

The cool star might also possibly be an FK Comae 
star (Bopp & Stencel 1981) in a wide binary with the 
ionizing star. If this is the case, intense Call H and 
K and MgTI emission is again predicted, as is direct 
evidence of rapid rotation. Unlike the RS CVn stars, 
no large radial velocity variations would be expected. 

An accurate systemic radial velocity will also pro- 
vide the necessary constraints on the space motion 
so that K 1-6 can be definitively classified as either 



a disk or halo PN. In the future, an accurate paral- 
lax will come from the 

GAl43 

astrometric satellite. 
This space mission, with its microarcsecond positions, 
proper motions and parallaxes, along with its ability to 
determine accurate radial velocities, promises to revo- 
lutionise our understanding of the Milky Way's stellar 
content, PNe included (see van Leeuwen 2007, for a 
fuller discussion). 

6 Conclusions 

We obtained narrow-band images of the little studied 
emission nebula, K 1-6 with the FTN. A full search 
using online tools uncovered a range of multiwave- 
length archival data on the nebula and its central star. 
The study was conducted in conjunction with Year 11 
high school students as part of an Australian Research 
Council Linkage Grant science education project, de- 
noted Space To Grow. As Percy et al. (2008) have 
emphasised, an important goal is to allow students to 
"learn science by doing science" with real data sets. 

Our FTN Ha (+ [Nil]) images show K 1-6 to be 
a faint elliptical nebula with dimensions 197" x 160" 
that shows an asymmetric morphology indicative of 
an interaction with the ISM. As detailed in Section 
4.2, the available body of evidence (Frew & Parker 
2010) suggests the nebula is best interpreted as an 
evolved, low-surface-brightness PN based on its limb- 
brightened bow-shocked morphology, centrally concen- 
trated [O III] emission, ionized shell mass, vertical dis- 
tance from the Galactic plane, and the observed bright- 
ness of its ionizing star. 

The ofi^-centred "central" star has the photometric 
colours of a G/K subgiant or giant star, while archival 
GALEX images show a very hot subdwarf or pre- white 
dwarf coincident in position with this cool star. The 
central star is unusual, being strongly variable in the 
optical with a period of 21.312 ± 0.008 days and an 
amplitude of 0.7 mag. It inhabits a position in period- 
amplitude space that appears to be unique amongst 
variable PN central stars (cf. De Marco et al. 2008). 

Based on the stellar colours and X-ray properties, 
we prefer a RS CVn classification for the cool star, but 
the amplitude is unusual for the class. The presence 
of a hot ionizing source in the system implies that the 
system may be ternary, as most RS CVn systems con- 
tain two cool stars. Alternatively, the star is a member 
of the V471 Tan subclass, which is a binary containing 
an evolved cool star and a hot white dwarf. A system 
containing a FK Com star and a wide hot companion is 
also a possibility. Assuming a distance of 1.0 kpc, the 
inferred visual magnitude of the ionizing star is con- 
sistent with known PN central stars (Frew & Parker 
2010) and it appears to be located on the white dwarf 
cooling track in the HR diagram. 

The high transverse velocity of the binary central 
star marks this as a kinematically old system. The cool 
star must have been spun up to rotate at a sufficient 
speed to show strong X-ray emission, which is normally 
a signature of stellar youth. The cool star has either 
been spun up by mass transfer in a close binary or 

^^http://www.rssd.esa.int/gaia/ 
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via wind accretion in a wider binary system. Hence, 
the companion should show signs of enrichment in s- 
process elements such as barium, that were derived 
from the AGB wind; the central star of the evolved 
PN Abell 70 (Miszalski et al. 2010b) is just such a 
case. However, if the cool component is an FK Com 
star formed by a binary merger, then no enrichment 
signature will be seen. 

We urge further observations of this interesting sys- 
tem, including deep spectroscopic follow-up of the neb- 
ula and in particular, medium resolution spectroscopy 
of its central binary (or ternary) star in order to re- 
move any ambiguity in its classification. 
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